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FOREWORD 

The  Defense  Communications  Engineering  Center  (DCEC)  Technical 
Notes  (TN's)  are  published  to  inform  interested  members  of  the  defense 
community  regarding  technical  activities  of  the  Center,  completed  and 
In  progress.  They  are  intended  to  stimulate  thinking  and  encourage 
information  exchange;  but  they  do  not  represent  an  approved  position 
or  policy  of  DCEC,  and  should  not  be  used  as  authoritative  guidance  for 
related  planning  and/or  further  action. 

Comments  or  technical  Inquiries  concerning  this  document  are  welcome, 
and  should  be  directed  to: 

Director 
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I.  INTRODUCTION 


-'-—^In  order  to  evaluate  the  potential  impact  of  facsimile-imagery 
service  upon  the  current  and  future  Automatic  Voice  Network  (AUTOVON)  and 
to  reduce  the  time  required  to  perform  worldwide  AUTOVON  Flash  Non-Blocking 
(FNB)  Studies,  the  capability  of  the  DCEC  Voice  Network  Model  has  been 
extended  from  a single  class  traffic  model  to  a two  class  traffic  model. 

As  in  the  AUTOVON  system,  one  class  of  traffic  can  preempt  the  other  class 
of  traffic.  Both  ruthless  and  friendly  preemption  disciplines  have  been 
incorporated  in  the  model  for  either  originating  office  or  spill  control. _ 
This  report  documents  the  approach  used  by  DCEC  to  extend  the  DCEC  Voice// 

kux.  u.j.i  j.  _ . a-.  . - i. rx:  - - i j-i  . • l. a.  j " ' ^ 


Network  Model  to  a two  traffic  class  model  with  preemption. 
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II.  BACKGROUND 


DCEC's  Voice  Network  Model  Cor  algorithm)  has  been  modified  to 
include  traffic  capable  of  preemption.  The  details  of  the  modification 
of  this  algorithm  have  been  documented  in  this  technical  note.  The 
algorithm  was  initially  developed  by  Katz  [1],  implemented  at  DCEC  by 
Fischer  and  Knepley  [2],  and  incoporated  into  an  overall  voice  network  design 
and  analysis  algorithm  (see  Fischer,  Garbin,  Harris,  and  Knepley  [3])., 

It  is  an  analytical  model  hased  on  queueing  theory  which  converges  to 
a steady  state  description  of  the  network. 

Since  this  algorithm  provides  the  basis  from  which  we  start,  it 
will  be  described  briefly.  For  a more  detailed  discussion,  see  refer- 
ences [1],  [2],  and  [3].  The  algorithm  assumes  that  the  network 
topology  is  known.  That  is,  the  number  and  locations  of  the  nodes  of 
the  network,  the  links  and  channel  cross  sections,  point  to  point 
traffic  requirements,  and  routing  philosophy  are  all  given.  It  also 
assumes: 

a.  Traffic  originates  at  the  nodes  of  the  network  as  a Poisson 
process.  For  the  preemptive  network  we  assume  that  both  classes  of 
traffic  originate  in  a Poisson  process. 

b.  The  length  of  time  required  to  serve  a call  can  be  described 
by  an  exponential  distribution.  For  the  preemptive  network  we  assume 
both  classes  of  traffic  have  the  same  exponentially  distributed  holding 
time. 

c.  The  network  is  in  steady  state, 

d.  The  mean  and  variance  of  the  offered  load  are  sufficient  to 
describe  the  behavior  of  a customer.  This  assumption  is  only  used  in 
the  ruthless  preemption  discipline. 

e.  The  behavior  of  calls  on  different  links  of  the  network  is 
statistically  Independent. 

f.  Congestion  in  the  network  is  only  caused  by  the  lack  of  avail- 
able channels  on  the  links. 

In  this  network  no  queueing  of  calls  is  allowed  and  calls  may  search 
several  paths  In  the  network  until  they  find  a free  one.  In  trying  to 
mathematically  describe  what  is  happening,  one  is  forced  to  consider  a 
network  of  queues  where  at  each  queue  no  waiting  is  allowed  and  calls 
may  have  several  paths  to  their  destination  from  the  same  origin. 
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Some  of  the  typical  outputs  that  are  available  are  shown  in  Figure 
1.  The  same  outputs  are  available  in  the  modified  algorithm  for  each 
class  of  traffic. 


OUTPUTS 


NETWORK  N0DE-T0-N0DE  LINK 


GRADE  OF  SERVICE  X 

LINK  UTILIZATION 
LINK  GRADE  OF  SERVICE 
CONGESTION  POINTS  X 


X 

X 

X 

X 


Figure  1.  Typical  Outputs  of  Voice  Performance  Algorithm 


Typical  numerical  and  computer  time  comparisons  with  tne  results  of  an  event 
by-event  simulation  are  given  in  reference  [3]. 

For  the  network  with  precedence  traffic  we  assume  there  are  two 
classes  of  calls  trying  to  use  the  network,  classes  1 and  2.  Class  1 
calls  have  priority  over  the  class  2 calls.  Both  classes  of  calls  use 
the  same  routing  tables.  Two  types  of  preemption  disciplines  are 
considered.  The  first  is  known  as  ruthless;  here,  a class  1 call 
finding  all  the  channels  on  a link  busy  preempts  a class  2 call  that 
occupies  a channel.  If  no  class  2 calls  are  present  (i.e.,  class  1 calls 
are  using  all  of  the  channels)  the  arriving  class  1 call  tries  to  find 
an  alternate  path  in  the  network.  Preempted  class  2 calls  are  assumed 
lost.  Class  2 calls  arriving  at  a link  try  an  alternate  path  if  all 
channels  in  the  link  are  busy. 


The  second  type  of  preemption  discipline  that  is  considered  is 
called  friendly;  here  a class  1 call  arriving  at  a link  with  all  the 
channels  busy  tries  to  find  an  alternate  path  in  the  network  to  the 
destination.  If  all  alternate  paths  are  occupied,  then  the  class  1 call  re- 
turns to  the  oriqinal  link  and  preempts  a class  2 call  that  occupies  a channel. 
Again  if  all  the  channels  are  occupied  by  class  1 calls,  the  call  is 
considered  blocked  on  that  link  and  tries  an  alternate  route.  At  the 
next  link  the  same  discipline  is  used.  For  the  class  2 calls,  the  same 
philosophy  and  call  disposition  occurs  as  in  the  ruthless  case. 
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The  basic  mathematical  problem  that  one  encounters  is  that  the 
output  process  from  a link  is  not  Poisson  and  is  not  simply  described. 

This  is  further  complicated  by  the  fact  that  there  are  three  types  of 
traffic  trying  to  use  any  given  link:  originating,  tandem,  and  overflow. 
Originating  traffic  is  that  traffic  which  has  just  entered  the  network; 
tandem  traffic  is  that  traffic  carried  previously  by  another  link  of  the 
network;  and  overflow  traffic  is  that  traffic  which  has  been  blocked  and 
is  trying  to  find  an  alternate  path  to  the  destination.  The  arrival 
processes  for  each  of  these  are  different,  and  are  not  Poisson  for  the 
tandem  and  overflow  (see  Cooper  [4]).  Thus,  each  type  of  call  sees  a 
different  blocking  probability  (see  Kuczura  [5] ) and  the  evaluation  of 
this  probability  is  not  straightforward  (Kuczura  [5]  and  Disney  [6]). 

Thus,  one  is  forced  to  consider  some  sort  of  procedure  which  approximates 
the  behavior  of  the  network.  We  have  adapted  an  iterative  procedure, 
first  introduced  by  Katz  [1],  that  uses  the  mean  and  variance  of  the 
offered  load  to  describe  the  behavior  of  a call  on  a link. 

The  algorithm  is  an  iterative  procedure  which  is  based  on  the 
assumption  that  the  traffic  flow  in  the  network  can  be  adequately  des- 
cribed by  the  mean  and  variance  of  its  distribution.  This  two-parameter 
characterization  accounts  for  the  non-Poisson  nature  of  the  traffic  which 
overflows  or  is  carried  by  a link. 

At  each  iteration,  a two-step  procedure  is  implemented  in  the  algor- 
ithm. The  first  step  is  a load  assignment  procedure.  During  this  step 
the  offered  traffic  is  distributed  over  the  network  (one  source-desti- 
nation pair  at  a time)  using  the  routing  table  and  the  link  blocking  

probabilities  calculated  at  the  last  iteration.  Initially,  these 
probabilities  are  set  to  zero.  By  considering  all  source-destination 
pairs  of  the  network,  the  offered  load  to  each  link  of  the  network  and  the 
total  traffic  lost  for  all  pairs  are  accumulated  during  the  load  assign- 
ment procedure. 

The  next  step  in  the  current  iteration  is  to  compute  new  estimates 
of  the  link  blocking  probabilities  and  the  mean  and  variance  of  the 
overflow  processes.  From  the  load  assignment  procedure  we  know  the  mean 
and  variance  of  the  offered  load  to  each  link  of  the  network.  Based  on 
the  relationship  between  these  parameters,  three  methods  are  used  to 
determine  blocking  probabilities  and  moments  of  overflow  processes.  If 
the  mean  and  variance  are  equal,  we  then  assume  the  offered  load  is 
being  generated  by  a Poisson  process  and  Erlang's  Loss  Formula  is  used 
to  give  an  expression  for  the  blocking  probability.  For  the  case  where 
the  mean  is  less  than  the  variance,  the  process  is  basically  one  for 
which  the  traffic  has  overflowed  another  link  (see  Cooper  [4])  and 
Wilkinson’s  Equivalent  Random  Technique  (see  [4])  is  used  to  obtain  the 
desired  results.  For  the  case  where  the  mean  is  greater  than  the 
variance,  the  Smoother  Than  Random  Technique  (see  Katz  [1])  is  used. 
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If  the  new  estimates  are  not  sufficiently  close  to  the  computed 
values  of  the  previous  iteration,  the  load  assignment  procedure  is 
repeated  using  the  new  values  of  link  blocking.  The  algorithm  converges 
when  the  link  blocking  probabilities  computed  from  a set  of  offered 
loads  are  the  same  as  those  used  to  generate  those  loads.  Typically, 
convergence  requires  from  7 to  10  iterations. 


We  have  modified  the  basic  voice  performance  algorithm  to  handle 
either  of  the  preemption  disciplines  mentioned  above.  The  current 
version  of  the  algorithm  assumes  the  same  discipline  is  used  throughout 
the  network.  In  this  technical  note  the  philosophy  of  the  required 
modifications  is  discussed  (section  III).  In  section  IV  some  numerical 
and  time  comparisons  with  the  results  of  an  event-by-event  simulation 
are  given.  The  conclusions  and  future  work  are  discussed  in  section  V. 
For  completeness,  an  appendix  with  three  subsections  is  included.  In 
the  first  subsection  of  the  appendix  the  mathematical  analysis  of  the 
link  performance  of  the  friendly  discipline  is  given.  Subsection  2 
gives  a mathematical  treatment  of  a link  rea liability/ performance 
analysis  required  in  the  Flash  Non  Blocking  Study.  The  last  subsection 
contains  a detailed  discussion  of  the  load  reduction  procedure  used 
in  the  basic  network  design  algorithm. 
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III.  MODIFICATIONS  TO  THE  BASIC  ALGORITHM 


From  section  II,  one  can  see  that  a vital  ingredient  of  the  basic 
performance  algorithm  is  the  capability  to  predict  link  performance. 

Thus,  for  each  of  the  preemption  disciplines,  ruthless  and  friendly,  one 
must  be  able  to  predict  the  link  performance  for  each  class  of  traffic. 
Once  the  required  equations  have  been  developed,  the  appropriate  changes 
to  the  basic  algorithm  have  to  be  made.  In  this  section,  we  discuss  the 
link  models  and  the  subsequent  changes  required. 

First,  let  us  consider  the  link  performance  of  the  ruthless  pre- 
emption discipline.  We  assume  that  the  arrival  process  of  each  class  to 
the  link  is  an  independent  Poisson  process,  and  that  the  link  has  s 
channels.  The  load  for  class  i calls  is  Pi  erlangs  ( i = 1 , 2 ) . As  in 
section  II,  the  service  holding  times  are  exponentially  distributed  with 
the  same  mean  for  each  class  of  calls.  Class  1 calls  have  precedence 
over  class  2 calls  in  the  sense  that  if  all  the  channels  are  busy  when  a 
class  1 call  arrives,  the  class  1 call  will  preempt  a class  2 call  on  the 
link.  If  there  are  only  class  1 calls  on  the  link,  the  arriving  class  1 
call  is  alternately  routed.  The  preempted  class  2 call  is  assumed  lost 
and  is  not  alternately  routed.  Only  arriving  class  2 calls  which  find 
all  channels  busy  are  alternately  routed.  We  are  interested  in  deter- 
mining the  steady  state  loss  probability  for  each  class,  denoted  by 
PL-f ,1*1,2.  It  should  be  pointed  out  that  one  minus  this  measure  of 
performance  is  the  proportion  of  class  i traffic  that  completes  its  call 
on  this  link.  The  quantity  PL i does  not  always  represent  the  blocking 
probability  for  class  i on  the  link.  The  blocking  probability  is  the 
probability  an  arriving  call  finds  all  the  channels  busy. 

The  mathematical  analysis  of  this  preemption  discipline  has  been 
considered  by  Burke  [7].  If  E(a,s)  is  Erlang's  loss  formula,  that  is 


E(a,s) 


then 

PL-,  = E(Pl,s),  (1) 


aS/ 


s! 


I an/ 


n=0 


n! 
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. : 


and 


p-| {E(p1+p2,s)-E(p-|  ,s)} 


PL 


2 


E(pi+p2,s)  + 


p2 


(2) 


Equation  (1)  is  easy  to  see,  since  the  class  1 calls  preempt  class  2 
calls  when  blocked;  the  only  traffic  class  1 calls  have  to  contend  with 
are  from  its  own  class.  Erlang's  Loss  Formula,  £(a,s),  represents  the 
probability  of  s busy  channels  with  an  offered  load  of  a erlangs.  The 
only  time  a class  1 call  does  not  receive  service  from  the  link  (and 
hence  is  lost)  is  when  all  the  channels  are  busy  with  class  1 calls. 
Equation  (1)  then  follows. 

Equation  (2)  requires  a little  more  explanation.  Class  2 calls  are 
denied  service  from  the  link  in  two  ways.  First,  all  the  channels  may 
be  busy  on  the  initial  attempt,  and,  second,  a class  2 call  may  be  pre- 
empted after  it  has  started  to  use  one  of  the  channels.  It  can  be  shown 
that  the  probability  that  all  s channels  are  busy,  regardless  of  the  class 
of  calls  using  the  link,  is  given  by  Erlang's  Loss  Formula  with  the 
offered  load  equal  to  the  sum  of  the  class  1 and  class  2 loads;  i.e., 
a=P!+P2.  This  factor  is  given  by  the  first  term  in  the  right-hand 
side  of  equation  (2). 

The  second  way  a class  2 call  does  not  receive  service  from  the 
link  is  by  preemption.  This  factor  is  given  by  the  second  term  in  the 
right-hand  side  of  equation  (2);  but  let  us  look  at  it  more  closely. 

The  quantity  E(pi+P2>s)-£(pi,s)  represents  the  probability  that  class 
1 calls  preempt  class  2 calls;  this  factor  multiplied  by  the  class  1 
load,  pi,  gives  the  class  2 load  that  is  preempted.  The  whole  quantity 
divided  by  the  offered  class  2 load,  P2,  is  the  portion  of  class  2 calls 
that  are  preempted  on  the  link.  This  quantity  is  the  second  term  in  the 
right-hand  side  of  equation  (2). 

From  earlier  discussions  on  the  effect  of  a ruthless  preemption 
discipline  on  class  1 calls,  one  can  see  that  there  is  no  change  in 
the  basic  algorithm  to  determine  the  network  performance  for  the  class 

1 calls.  In  fact,  one  may  use  the  two-moment  version  of  the  alaorithm 

to  describe  the  behavior  for  class  1 calls.  For  all  the  other  cases  that 
we  consider,  only  a one-moment  characterization  of  the  offered  traffic 
is  used. 

For  class  2 calls  changes  are  required.  First,  equation  (2)  is 
used  at  the  end  of  the  Iteration  to  update  the  link  measures  of 
performance,  based  on  the  current  estimate  of  the  offered  class  1 and 

2 loads.  During  the  load  assignment  procedure  for  class  2 traffic, 
there  are  some  additional  changes  to  be  made,  pictorially  represented  in 
Figure  2. 
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Figure  2.  Load  Assignment  Procedure  for  Class  2 Calls 


Let  p^,i=l ,2,  be  the  estimated  class  i load  on  link  (k,£)  from  the 

j iteration,  and  su  » the  number  of  channels  on  link  (k,£).  In  the 

st 

(j+1)  iteration,  suppose  a class  2 call  wants  to  go  from  node  k to 
node  m,  and  the  routing  dictates  that  blocked  calls  on  l-’nk  (k,*) 
overflow  to  link  (k,n).  Class  2 calls  would  overflow  to  link  (k,n)  with 
probability  E(pi+p>ltsk*),  while  class  2 calls  would  be  carried  by  link 
(k,i)  with  probability  1-PL2H,pl,sk  z ),  where  PL2(>i,p^,sk  4 ) is 
equation  (2)  with  parameters  p^,p^,  and  sk  % . That  is,  class  2 calls 
are  overflowed  via  the  blocking  (all  channels  busy)  formula  and  tandemed 
with  the  probability  it  is  carried  by  link  (k,z).  We  note  from  equation 
(2)  that  PL2 (p i ,P2*sia)  is  greater  than  or  equal  to  E(PJ+PJ,sk  A). 

For  the  friendly  preemption  discipline,  more  problems  arise.  Since 
a blocked  dies  1 call  tries  to  find  another  path  in  the  network  before 
searching  the  link  in  a ruthless  manner,  the  link  performance  analysis 
is  difficult.  In  order  to  make  the  mathematical  analysis  tractable,  we 
formulate  the  problem  as  show  in  Figure  3. 


Figure  3.  Link  System  Formulation  for  Friendly  Search 

i 
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Let  Pn- i = 1 ,2  be  the  offered  class  i load  to  a link  (k,2)  with  SkfA 
channels,  and  the  statistical  independence  assumptions  be  the  same  as 
in  the  ruthless  case.  An  arriving  class  1 call  that  finds  all  the 
channels  busy  tries  an  alternate  path  in  the  network.  If  all  the  paths 
are  occupied,  it  returns  and  searches  link  (k,t)  in  a ruthless  manner. 
We  assume  that  the  probability  of  returning  is  6.  The  complete  mathe- 
matical analysis  of  this  system  is  given  in  the  appendix. 

We  are  interested  in  determining  the  loss  probability,  PLi,  for 
each  class  of  calls  If  Q]  represents  the  number  of  class  1 calls  on 
the  link  in  steady  state,  then 

PL  "I  = (1”6)E(pi+P2*sk^)+SPr{Qi=$k^} 


and 


P,  - rf  , - s,  ple{E(pl+p2’Sk,Jl)-Pr{Ql=sk,t)}} 

PL2  - E(Pl+p2,sK>Jl)+  • 


Let  us  first  consider  the  loss  probability  for  the  class  1 calls,  PL]. 
The  quantity  E(p1+p2,sk  ) represents  the  probability  all  the  channels 
are  busy  and  hence  theKprobability  an  arriving  class  1 call  tries  to 
find  an  alternate  path  in  the  network.  Since  1-6  is  the  probability 
the  call  does  not  return  to  the  link,  the  first  part  on  the  right- 
hand  side  of  equation  (3)  is  the  portion  of  arriving  class  1 calls 
that  find  all  the  channels  busy  and  a free  alternate  path.  These  calls 
are  considered  lost  to  link  (k,t). 


There  is  another  way  calls  can  be  lost  to  the  link.  The  prob- 
ability PrlQ^s^  » }is  the  probability  the  channels  are  all  occupied  by 
class  1 calls.  If  we  multiply  this  quantity  by  the  probability  all  the 
alternate  paths  are  occupied,  6,  the  product  represents  the  portion 
of  class  1 calls  that  find  all  the  channels  on  the  link  occupied  by 
class  1 calls  and  all  the  alternate  paths  busy.  The  quantity  is  given 
by  the  second  part  in  equation  (3). 


The  development  of  the  loss  probability,  PL2.  for  class  2 calls  is 
similar  to  the  case  of  ruthless  preemption.  The  only  quantity  that  needs 
some  explanation  is  the  numerator  of  the  second  factor  in  equation  (4). 
The  quantity  E(<"1+P2»sk  t)-Pr{Qi*S|<t4}  is  the  probability  that 
class  1 calls  preempt  class  2 calls!  This  quantity  multiplied  by  the 
amount  of  Class  1 traffic  capable  of  preempting  class  2 traffic,  6plt 
is  the  amount  of  preempted  class  2 traffic.  Dividing  this  quantity  by 
the  class  2 traffic  offered  to  the  link  gives  the  probability  that  a 
class  2 call  will  be  preempted.  The  quantities  Pj'and  p2  are  used  as 
they  were  in  the  ruthless  case. 

Equations  (3)  and  (4)  are  used  In  the  algorithm  in  exactly  the 
same  way  as  equations  (1)  and  (2)  were  in  the  ruthless  preemption  dis- 
cipline. So  changes  to  the  basic  algorithm  follow  along  the  same  line  as 


discussed  for  the  ruthless  discipline.  There  still  remain  two  problems 
r before  equations  (3)  and  (4)  can  be  implemented.  We  have  not  given 

an  easily  computable  expression  for  Pr{Qi=Skti}.  How  does  one 
determine  0 in  the  algorithm?  as  we  mentioned  earlier,  a complete  mathe- 
matical analysis  of  the  link  performance  is  given  in  the  appendix.  From 
V that,  one  can  determine  Pr{Q-]=s}  (note  that  we  have  dropped  the  sub- 

scripts indicating  the  link  under  consideration)  from  the  solution  of  s-1 
linear  equations.  This  type  of  solution  would  not  be  feasible  for  imple- 
mentation in  our  algorithm  because  it  would  require  too  much  computer 
run  time. 

One  is  then  forced  to  consider  some  sort  of  approximation  for 
PrlQj  = s}.  The  details  of  the  approximation  are  given  in  the  appendix; 
but  the  basics  are  as  follows.  We  have  examined  some  numerical  data 
where  we  have  numerically  solved  for  the  correct  PKQ-|=s  }.  This 
investigation  indicated  that  Pr(Q^=s}  had  a quadratic  from  in  6.  Let 
us  denote  the  approximation  to  Pr{  Q-j=s } by  f(pi,P2>s,3)  and  thus 

f ( Pi » » s > S ) =^o+^l B+C26  , (5) 

where  the  constant  Cq,  C]  and  C2  have  to  be  determined. 

When  3=0  the  blocked  class  1 calls  do  not  return  and  the  behavior  on  the 
link  is  as  if  both  classes  of  calls  were  equally  contending  for  the 
channels.  From  Cooper  [4],  one  can  show 

s 

Pl/  1 

Pr(Q,-s)  • — . (6) 

j0<OlVn/n! 

Similarly,  if  0=1  then  the  class  1 calls  appears  as  if  they  were 
operating  in  the  ruthless  mode  and 

Pr{Q1=s}  = E(pvs).  (7) 

Equations  (6)  and  (7)  provide  two  of  three  equations  necessary  for 
determining  the  three  unknowns  (Cq,  C],  C2);  the  other  equation  is 
found  by  examining  numerical  data  and  relating  the  parameters  to  the 
boundary  condition  0=1: 


(p-j+p2)f(pi  ,p2»s.0.6)%pif(pi  ,p2fs,l ) (8) 

where  X means  approximate. 


r 
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Equations  (6),  (7),  and  (8)  provide  three  equations  in  Co,  Cl,  and  Cz . 
The  solution  of  these  equations  and  the  comparison  of  the  approximation 
to  the  exact  solution  are  given  in  the  appendix. 


The  other  problem  is  the  determination  of  8 for  each  link  of  the 
network.  Remembering  that  8 represents  the  probability  of  an  arriving 
class  1 call  returning  to  the  link  after  searching  all  the  alternate 
paths  in  a friendly  mode,  we  essentially  set  8 for  each  link  to 

8 = keAE(P1(k)+P2(k),Sk)  , (9) 

where  A is  the  set  of  alternate  links  for  the  link  under  consideration, 

Pi(k)  the  class  i load  on  the  kth  link,  and  sk  the  number  of  channels  in 
the  link.  In  reality,  8 should  be  d i recti onali zed;  that  is,  calls 
going  from  node  k to  node  i have  a different  6 than  calls  from  node  i to 
k.  This  fact  has  not  been  incorporated  into  the  current  version  of  the 
algorithm.  At  the  end  of  each  iteration  a new  6 is  computed  for  each 
link. 

In  the  next  section,  the  current  version  of  the  voice  preemption 
model  is  compared  with  the  results  of  an  event  by-event  simultation. 


IV.  NETWORK  MODEL  PERFORMANCE 


The  performance  of  the  DCEC  Voice  Network  Model  with  Preemption  was 
determined  for  both  the  ruthless  and  friendly  preemption  disciplines. 

For  each  preemption  discipline,  five  event-by-event  computer  simulation 
runs  (each  using  a different  random  seed)  were  performed  for  three 
different  network  configurations  representative  of  the  Pacific,  CONUS, 
and  European  AUTOVON.  Next,  using  the  traffic  data  from  the  simulation 
runs  as  input  to  the  DCEC  Voice  Network  Model  with  Preemption,  the 
network  and  switch-to-switch  grades  of  service  (GOS's)  were  estimated  by 
the  model  for  each  configuration.  The  95%  confidence  intervals  for  the 
network  and  switch-to-switch  GOS's  obtained  from  the  computer  simulation 
runs  were  computed  using  Student's  t distribution.  The  GOS's  estimated 
by  the  model  were  then  compared  to  the  average  GOS's  for  the  five 
computer  simulation  runs  and  the  number  of  estimates  from  the  model 
which  fell  within  the  confidence  interval  for  the  computer  simulation 
runs  was  determined. 

Figure  4 shows  the  results  of  the  analysis  of  the  performance  of 
the  DCEC  Voice  Network  Model  with  Preemption  for  the  ruthless  preemption 
discipline.  For  each  of  the  three  configurations  over  three  fourths  of 
the  switch-to-switch  GOS's  fell  within  the  95%  confidence  interval 
derived  from  the  results  of  the  computer  simulation  runs  of  the  perfor- 
mance of  each  network  configuration. 

Since  the  network  GOS  for  first  class  traffic  was  less  than  or 
equal  to  0.05  for  the  CONUS  and  European  configurations,  a dual  moment 
model  for  class  1 traffic  was  used  in  order  to  maintain  the  accuracy  of 
the  model.  (Below  a network  GOS  of  0.05  a dual  moment  model  is  usually 
necessary  in  order  to  maintain  accuracy.) 

Figure  5 shows  the  results  of  the  analysis  of  the  performance  of 
the  DCEC  Voice  Network  Model  with  Preemption  for  the  friendly  preemption 
discipline.  With  the  exception  of  the  CONUS  configuration,  over  three 
fifths  of  the  switch-to-switch  GOS's  fell  within  the  95%  confidence 
interval  derived  from  the  results  of  the  computer  simulation  runs  of  the 
performance  of  each  configuration.  For  the  Pacific  network  configuration, 
dual  moment  modeling  was  not  necessary  since  the  network  GOS  for  each 
class  of  precedence  traffic  exceeded  0.05.  Using  a single  moment  model 
for  CONUS,  only  one  third  of  the  number  of  switch-to-switch  GOS's  fell 
within  the  95%  confidence  Interval  for  class  1 traffic  (the  overall 
network  GOS  equaled  .04).  No  dual  moment  model  as  yet  has  been  imple- 
mented for  the  friendly  preemption  discipline.  Although  the  European 
network  also  experienced  a network  GOS  less  than  0.05  for  class  1 traffic 
the  percentage  of  switch-to-switch  pairs  whose  GOS  fell  within  the  95% 
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confidence  interval  was  oyer  three  fourths  (82%),  a considerable  improvement 
over  the  CONUS  configuration  (32%).  This  improvement  was  accomplished  by 
incorporating  a 5%  increase  of  the  non-terminating  tandem  class  1 traffic  carried 
on  each  interswitch  trunk  in  the  model  in  order  to  account  for  peaicedness  of 
any  overflow  traffic. 

_Figure  6_shows._a__cj3mpa_ri son  Qf_the  central  processing  unit  (CPU)  time 
of  an  IBM  370/155  computer  used  during  simulation  and  computer  runs  nf 
the  Voice  Network  Model  with  Preemption. ; 


AREA 

SIMULATION 

KVMggglgjmB 

Pacific 
(6  Switches) 

23  min. 

8 sec. 

CONUS 

(10  Switches) 

35  min. 

37  sec. 

Figure  6.  Comparison  of  IBM  370/155  CPU  Time  Used  by  Event-by-Event 
Simulator  and  the  Voice  Network  Model  with  Preemption 


In  order  to  minimize  computer  turn-around  time,  five  short  runs  were  made 
for  each  simulation.  A total  of  216  thousand  calls  v/ere  processed  for 
the  Pacific  and  338  thousand  for  the  CONUS.  A new  random  seed  was 
selected  for  each  run  in  order  to  eliminate  any  possible  bias  in  the 
generation  of  traffic. 

The  CPU  time  to  determine  the  performance  of  the  full  CONUS  AUTOVON 
configuration  (60  switches)  using  the  unmodified  Voice  Network  Model  for 
a single  class  of  traffic  (5000  erlangs  of  ROUTINE  traffic)  is  about  10 

to  12  minutes.  Since  the  Voice  Network  Model  with  Preemption  may  be 

viewed  as  two  unmodified  models  running  in  parallel,  one  for  class  1 
traffic  and  the  other  for  class  2 traffic,  the  CPU  time  for  the  Voice 
Network  Model  with  Preemption  for  the  full  CONUS  AUTOVON  configuration 
should  be  no  more  than  twice  that  of  the  unmodified  model  or  20  to  24 
minutes.  If  the  precedence  traffic  is  much  less  than  that  of  the  ROUTINE 
traffic  then  the  load  assignment  procedure  for  the  precedence  traffic 
will  take  considerably  less  time.  There  are  about  200  to  300  erlangs  of 
precedence  traffic  in  the  CONUS  AUTOVON  during  the  busy  hour.  Thus  it 

may  be  estimated  that  the  Voice  Network  Model  with  Preemption  would 

use  about  18  CPU  minutes  to  estimate  performance  of  the  full  CONUS 
AUTOVON  configuration  during  the  busy  hour  (12  minutes  for  ROUTINE 
traffic  and  6 minutes  for  precedence  traffic).  Alternatively,  it  would 
not  be  unreasonable  to  estimate  that  the  corresponding  simulation  would 
require  about  three  CPU  hours.  Hence,  one  can  see  the  substantial 
savings  attainable  with  this  analytic  model. 
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V.  CONCLUSION  AND  FUTURE  WORK 


1 

A Voice  Network  Model  with  Preemption  has  been  developed  by  DCEC  which 
will  be  used  to  perform  future  engineering  studies  of  the  AUTOVON  involving 
preemption.  These  studies  could  only  have  been  performed  previously  by  J 

lengthy  computer  simulation.  Presently,  the  model  will  be  used  to  assess 
the  potential  performance  and  economic  impact  of  facsimile-imagery  service 
upon  the  present  and  future  AUTOVON.  It  will  also  be  used  to  reduce  the 
time  required  for  performing  FNB  studies.  All  AUTOVON  interswitch  trunking, 
which  must  meet  the  FNB  criteria,  will  initally  be  resized  using  the  model  ! 

in  place  of  event-by-event  computer  simulation. 

In  order  to  increase  the  usefulness  of  the  Voice  Network  Model  with 
Preemption,  the  following  extensions  to  the  model  are  planned: 

« 

(1)  Extend  the  capability  of  the  model  to  include  both  ruthless  and 
friendly  preemption  disciplines  in  the  same  network. 

(2)  Investigate  and  extend  the  model  to  include  different  holding 
times  for  each  class  of  traffic. 

(3)  Examine  the  precision  of  the  model  for  networks  similar  to  the 
ones  generated  in  the  FNB  study. 

(4)  Extend  the  capability  of  the  model  to  include  link  reliability. 

(5)  Develop  a network  design  program  for  two  classes  of  traffic  which 
woul  (^complement  the  Voice  Network  Model  with  Preemption. 

Both  ruthless  and  friendly  preemption  disciplines  exist  in  today's 
European  AUTOVON;  friendly  preemption  exists  intra-Europe  and  ruthless 
preemption  exists  at  the  gateway  switches  between  Europe  and  CONUS. 

Extending  the  capability  of  the  model  to  include  both  ruthless  and 
friendly  preemption  disciplines  in  the  same  network  model  will  enable 
DCEC  to  more  accurately  assess  the  performance  of  the  European  AUTOVON. 

In  general,  the  holding  times  of  each  class  of  traffic  may  not  have  equal 
distributions  as  assumed.  Thus,  the  effect  of  unequal  holding  times  should 
be  ascertained  and  incorporated  in  the  model. 

' 

Based  upon  analyses  performed  to  date,  the  precision  of  the  Voice 
Network  Model  with  Preemption  Is  adequate  for  estimating  network  grades 
of  service  of  not  less  than  one  or  two  blocked  calls  per  100  calls 
offered.  One  of  the  FNB  criteria  for  AUTOVON  is  that  no  more  than  one 
FLASH  call  per  1000  FLASH  calls  offered  can  be  blocked.  Thus,  in  order 
to  determine  whether  a network  meets  this  criteria  using  the  Voice 
Network  Model  with  Preemption,  the  precision  of  the  model  will  be 
studied  and  extended,  if  possible. 
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Appendix  A. 2,  Link  Performance  for  Unreliable  Channels,  outlines  the 
theory  necessary  to  extend  the  capability  of  the  model  to  include  link 
reliability.  A link  reliability  of  99%  is  used  in  the  computer  FNB 
simulation  runs.  The  Voice  Network  Model  with  Preemption,  with  the  link 
reliability  feature,  may  be  useful  in  performing  FNB  studies  as  well  as 
other  network  studies  involving  network  reliability. 

Finally,  we  would  like  to  incorporate  the  model  into  the  overall 
design  and  analysis  algorithm  currently  being  used  at  DCEC.  This  will 
extend  DCEC's  capability  to  design  as  well  as  to  analyze  communication 
networks  capable  of  preemption. 
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APPENDIX 


There  are  three  subsections  in  this  appendix.  The  first  deals 
with  the  mathematical  analysis  and  development  of  the  link  performance 
approximation  for  the  friendly  preemption  discipline.  In  the  second, 
a method  for  computing  link  performance  when  some  of  the  channels  in 
the  link  are  not  completely  reliable  is  discussed.  Finally,  an 
example  of  the  load  reduction  process  used  in  the  algorithm  is  given 
in  the  third  subsection. 

A. 1 LINK  PERFORMANCE  OF  FRIENDLY  PREEMPTION  DISCIPLINE 

Under  this  discipline,  an  arriving  high  precedence  (class  1)  call 
that  is  blocked  on  a link  first  tries  to  find  an  alternative  path  to 
the  node  at  the  other  end  of  the  link.  If  one  exists,  the  call  takes 
that  path;  if  not,  the  call  returns  to  the  link  and  preempts  a 
randomly  selected  lower  precedence  (class  2)  call.  If  no  class  2 calls 
are  present,  the  class  1 call  is  alternately  routed. 

It  is  assumed  there  are  s channels  in  the  link;  calls  from  each  of 
the  two  classes  arrive  at  the  link  in  accordant  with  two  independent 
Poisson  processes  with  parameters  Xi  and  A2.  respectively.  The 
holding  time  for  each  class  of  calls  is  assumed  to  be  the  same  and 
exponential ly  distributed  with  mean  u.  Arriving  class  1 calls  finding 
all  the  channels  busy  leave  the  link  and  try  to  find  another  path  in 
the  network.  We  assume  there  is  a probability,  e,  tha^  they  will  return. 
If  the  returning  class  1 call  finds  all  the  channels  occupied  by  class 
1 calls,  it  is  lost  to  the  link.  No  queueing  is  allowed  for  either 
class.  In  the  steady  state  condition,  we  are  interested  in  determining 
the  loss  probabilities  for  class  1 and  2 calls. 

In  steady  state,  let  , i =1 , 2 , be  the  random  variable  representing 
the  number  of  class  i calls  on  the  channels,  and  Q=Qi+Q2  the  total 
number  of  calls  on  the  channels.  If  p i i /v* . i=1.2,  and  pi  ,i=Pr{Ql=i , 

Q2= j > • then  the  steady  state  equations  for  P-jj  are  for  i=0j s-1 

and  j=0,l,2 s-l-i: 


(pi+P2+i+j)Pi  ,j=0ipi-i  ,j+(i+1  )pi+i  ,j+(j+1  )pi,j+l+p2pi  ,M  . 


(A.  1 ) 


and 


(piB+s)pijS.i=PiPi.i,s.i+p2Pi,s-i-l+Pl6Pi-l,s-i+l; 


sPs,O=PlPs-l,O+pi0Ps-l,l  , 


(A. 2) 
(A. 3) 


A-l 


where 


=Ps+l ,j=Pi 


For  the  development  of  birth  and  death  equations,  see  [4]. 

For  k=0,1,...,s,  let  R|<=Pr{Q=k};  then  from  equation  (A.1),  it  can 
be  seen  that  Rj<  must  satisfy  the  following  set  of  equations: 


(p1+p2+k)Rk=(p1+p2)R|<_1+(k+l  )Rk+1 , 

where  R_i=Rs+|sO.  The  solution  to  these  equations  is,  for  k=0,l,..,s 


Rk=Pr{Q=k)  = 


(p i+p2) 
k! 


s (Pi+P2)d 
" j! 


j=0 


(A. 4) 


(A. 5) 


From  equation  (A. 5),  it  can  be  seen  that  the  steady  state  probability 
that  the  number  of  calls  on  the  channel  is  k does  not  depend  on  0,  and 
hence,  the  total  carried  load,  CL,  is  also  invariant  under  B and  is  given 
by 

s 

CL  = Z kR.  = (p i+p2) (1-R  ) . (A. 6) 

k=Q  K s . 


By  rewriting  equation  (A.l),  one  sees  that  the  solution  of  equations 
(A. 1 )-(A.3)  depends  on  s+1  unknowns  Pq  q,Pq  ^.....Pq  From  Rq-Pq  q 
and  equation  (A. 3)  the  solution  can  be  further  reduced  to  the  s-1  unknowns 

Pq  rP0  2 Pq  s_i*  0ne  may  then  use  equation  (A. 2)  to  obtain  s-1 

equations  in  these  unknowns.  These  equations  were  solved  directly  by  an 
iterative  method  in  order  to  develop  an  approximation  for  the  probability 
of  lost  class  1 calls. 

For  each  class  of  calls,  it  is  desired  to  determine  the  probability 
that  a call  from  that  class  will  be  lost  on  the  link.  The  loss  proba- 
bility, PL-|,  for  class  1 calls  is 

PL-j  = (l-s)Pr{Q=s}+BPr{Q1=s}.  (A. 7) 

where,  again, 


q=q1+q2. 


A-2 


The  first  part  on  the  right  hand  side  of  this  equation  represents  the 
situation  where  the  class  1 call  is  blocked  on  its  friendly  attempt  to 
use  the  channels  and  does  not  return  to  preempt.  The  second  part 
represents  the  preemption  situation.  For  this  case,  class  1 calls  are 
lost  on  the  link  when  they  return  and  find  all  the  channels  occupied  by 
class  1 calls. 

For  class  2,  the  loss  probability,  PL2>  is  given  by 

Pi0[Pr{Q=s}-Pr{Q,=s}] 

PL,=Pr{Q=s}+  1 • (A. 8) 

2 Pz 

The  first  factor  in  equation  (A. 8)  represents  the  situation  where  all 
the  channels  are  busy.  The  second  is  the  proportion  of  class  2 calls 
that  are  preempted.  In  this  expression,  the  numerator  is  the  number 
of  erlangs  of  class  1 calls  that  preempt  the  class  2 calls.  This 
number  must  equal  the  number  of  class  2 calls  that  are  preempted. 

Therefore,  the  proportion  of  class  2 calls  preempted  may  be  found  by 
dividing  this  amount  by  the  amount  of  class  2 offered  traffic,  p2;  thus 
the  second  factor  in  equation  (A. 8)  is  determined.  In  equations  (A. 7) 
and  (A. 3)  the  probability,  Pr{Q=s),  is  given  by  equation  (A. 5)  with  k=s; 
i.e.,  R$.  The  other  probability,  Pr{Q^=s>,  is  given  by  ps>g* 

Several  special  cases  are  of  interest.  First,  the  situation  where 
0=0;  i.e.,  class  1 and  2 calls  equally  contend  for  use  of  the  channels. 

For  this  case,  one  can  show  (see  [4])  tliat 

piS  (A. 9) 

P = A? — — 

Z (p  1"*"P2  ) ' -i  I 

j=0 

The  second  special  case  is  where  3=1.  This  situation,  is  equivalent  to 
the  protocol  where  an  arriving  class  1 call,  which  finds  all  channels  busy, 
immediately  preempts  a class  2 call  using  a channel.  This  protocol 
has  been  analyzed  in  ref.  [7J.  One  finds  that 


s 

pi/ 


s! 


s.O  s 
j*0 


2 PiJ/,t 


(A. 10) 


These  two  special  cases  serve  as  bounds  on  P n,  and  hence  on  PL,  and 
PL2.  s,u  1 


Another  special  case  is  where  s*l , From  equation  (A. 5),  one  can 
derive 

p . 1 

0,0  1+p 1+P2  * 

from  equation  (A.  3 ) 

p °2  - - ■ 

Ofl  (1+Pi8)(1+Oi+P2) 


and  finally,  using  equation  (A. 5), 

p1+p16(pi+p2) 

P1.0  (1+p1+p2)(l+PlB) 


From  these  results,  the  loss  probabilities  may  be  found: 

p1(l+&(pi+P2))+P20"8)) 

PL  = , 

( 1 +p l+p2 ) ( 1 1 8 ) 


and 


P1+P2+P1B(1+Pi+P2) 
(1+p 1+p  2 ) ( 1 +p le ) 


The  mathematical  analysis  presented  so  far  requires  one  to  solve  a 
set  of  linear  equations  to  determine  Ps,0*  The  use  of  a model  which 
would  solve  these  linear  equations  would  result  in  excessive  computer 
run  times.  A simple  approximation  has  been  developed  which  eliminates 
the  need  to  solve  the  set  of  linear  equations  and  in  turn  significantly 
reduces  the  computer  run  time. 

From  equation  (A. 7)  and  (A. 8),  one  sees  that  the  only  real  problem 
in  computing  PL  1 and  PL2  is  in  determining  Pr{Qi*s}S!Ps,0-  Thus,  one 
needs  to  approximate  this  probability  in  order  to  implement  these  results 
in  the  network  algorithm.  From  the  numerical  cases  that  were  examined, 
with  fixed  values  of  p-j,  P2,  and  s,  PS(o  appears  to  be  a quadratic 
function  In  0;  thus,  Ps,0»  can  be  approximated  by  the  function 
*(pl,p2,s,e),  which  is  given  by: 


where  the  constants  C_-(pi,P2,s),  i=0,l,2,  are  to  be  determined.  Two 
equations  in  these  unknowns  can  be  immediately  found.  As  was  pointed 
out  earlier,  when  6*0 


s,0  “ s 


Z (pi+p2)J/., 


and  when  8=1 


s.O  " s 


2 pi, 
j=0  ' j ! 

From  equation  (A. 11)  this  implies  that 


f(pi.P2.s,e)=C0(pi,p2,s)  = — — — 

1 ( P l+P2  )^/  •! 

j=0  J 

and  similarly 


Cq(pi»P2«s)+Ci (pi ,P2»s)+C-(p1,p2,s) 


PI5/,, 


* 01 

1=0  J ’ 


One  more  equation  is  still  required  for  determination  of  the 
constants  C-j  (pi . ^.s)  1*0, 1,2.  Examination  of  the  numerical 
cases  showed  that  at  8 *.6, 

f(Pl»P2»S,.6)  £ Pl  f(p1,p2,S,l) 

P1+P2 

where  £ signifies  approximately.  Using  equations  (A. 11)  and  (A. 14) 
one  finds: 

Pif(pi,P2.s,l>!(pi+p2)[C0(p1,p2,s)+.6C1(p1,p2,s)+.36C2(p1,p2,s)]. 


.V  "i.  '•  rf  * 


(A. 12) 


(A. 13) 


(A. 14) 


The  solution  for  C.fpj.Pg.s)  is  straightforward  and  is  given  by: 


Cq(pi»P2»s) 


S I 

Pi  /SJ 


5 i 

l (Pl+p 2)  /,| 
j=0  J ' 


C2(pi>P2»s)  = J Cq(p i»P2.s)  + 


. 3p2-.  2p ! 

•12(p1+p2) 


E(p 1 »s) , 


(A. 15) 


(A. 16) 


and 


Ci(pi»P2»s)  - E ( p 1 , s ) - Cq(pj,p2,s)  - C2(pi.P2>s) » (A. 17) 

where  E(pi.s)  is  Erlangs  Loss  Formula;  i.e.,  the  right  hand  side  of 
equation  (A.13).  For  the  cases  where  0=0  and  0=1,  one  sees  that  the 
approximation  is  exact. 

In  tables  A.l  through  A. 6 some  numerical  comparisons  of  the  exact 
results  to  the  approximations  are  given.  For  each  of  three  cases, 
s=l 0. 1 5 and  20,  the  effect  that  different  values  of  0 have  on  light 
and  heavy  loadings  were  examined.  For  a fixed  number  of  channels  and 
fixed  total  offered  load,  three  cases  were  considered.  The  first  case 
assumed  that  the  class  1 load  was  one-third  of  the  class  2 load,  the 
second  assumed  that  the  two  loads  were  equal,  and  the  last  that  the 
class  1 load  was  three  times  the  class  2 load.  For  each  of  these 
cases  and  different  values  of  0,  comparisons  of  the  exact  (E)  and 
approximate  (A)  results  are  given  for  P$  q,  PL^,  and  PL2. 

For  all  the  cases  presented  in  the  tables,  one  sees  that  the  approxi- 
mation estimates  PLi  and  PL2  very  well.  When  it  was  low  for  one  of 
the  loss  probabilities,  it  was  high  for  the  other.  Of  course,  this  was 
to  be  expected  since  the  total  carried  load  is  constant  and  the 
approximation  is  consistent  in  this  regard.  Another  comment  worthy  of 
mention  is  the  great  effect  0 has  on  PL]  and  PL2.  For  example, 
consider  Table  A. 2 with  pi=3  and  P2=9;  when  0*.2  PL^*.2415,  but  PLi= 

.0607  with  6*. 8.  It  can  be  seen  that  the  loss  probability  is  reduced 
by  a factor  of  four  when  0 is  increased  by  a factor  of  four. 

These  tables  certainly  demonstrate  the  effect  that  this  preemption 
discipline  has  on  the  performance  of  the  link.  It  tends  to  let  more 
class  2 calls  use  the  link  and  forces  the  class  1 calls  to  find  other 


f 


‘ 


I 

TABLE  A.  1 . LOW  LOAD  (p.+p,»6)  FOR  s*10 

1 * 


01=1.5,  p,=4.5 


TABLE  A. 2.  HIGH  LOAD  (pj+Pg-^)  FOR  s=10 


o,*3,  p,*9 


TABLE  A.  3.  LOW  LOAD  (t^+0^12)  FOR  S-15 

-9 


.0000 

.0000 


.0972 


.0000 

.0000 


.0000 

.0000 


.0000 


.0686 


.1028 


.0001 

.0002 


.0002 

.0004 


TABLE  A. 6.  HIGH  LOAD  Co-,  +p2»36)  FOR  s*20 


c-,=9,  p2=27 


paths  in  the  network.  This  fact  is  brought  out  in  Table  A. 2 for  the 
case  where  <M=9  and  P2=3.  In  the  ruthless  discipline  the  class  1 calls 
are  allowed  to  preempt  the  class  2 calls  immediately;  the  loss  proba- 
bility for  class  2 would  be  greater  than  .7.  But  under  the  friendly 
discipline,  if  B=.2,  then  Pl_2=.44.  Thus  more  class  2 calls  are  able  to 
use  the  link. 

A. 2 LINK  PERFORMANCE  FOR  UNRELIABLE  CHANNELS 

As  part  of  the  Flash  Non-Blocking  Study,  the  capability  is  required 
to  determine  the  link  performance  in  terms  of  loss  probability  when  the 
channels  on  the  link  fail.  Currently,  this  situation  is  handled  for 
the  Flash  Non-Blocking  Study  by  simulating  this  process  in  a Monte  Carlo 
manner.  That  is,  suppose  there  are  s channels  on  the  link  and  for  the 
time  period  under  consideration,  let  p be  the  probability  that  a channel  on 
the  link  is  working.  For  the  Flash  Non-Blocking  Study  and  for  each  channel  in 
a Monte  Carlo  run,  a uniform  random  number  is  drawn  and  compared  with  p in 
order  to  determine  if  the  channel  is  available.  Then,  the  number  of  channels 
available  on  each  link  is  used  in  an  event-by-event  simulation  of  the  system. 
After  the  simulation  is  in  a steady  state,  statistics  are  taken  and  the  process 
is  repeated.  That  is,  s more  uniformly  distributed  random  numbers  are  drawn 
and  compared  with  the  p's.  This  method  is  continued  until  some  statistical 
confidence  is  reached. 

In  this  appendix,  it  is  shown  how  to  solve  the  problem  mathematically. 
Again  let  s be  the  number  of  channels  and  p the  offered  load  in  erlangs, 
generated  by  a Poisson  process.  Assume  there  is  only  one  class  of 
customers  for  this  comparison.  If  p is  the  probability  a channel  is 
available  over  the  time  interval  of  interest,  then  the  probability  that 
there  are  n channels  available  is  given  by  a binomial  distribution.  If 
there  are  n channels  available,  then  Erlang's  Loss  Formula  E(p,n)  is  the 
equation  for  the  loss  probability.  Using  conditional  probability 
statements  one  arrives  at  the  loss  probability,  PL,  for  the  system  as: 


PL  = Z (p)pn(l -p)S_nE(p,n) . (A. 12) 

n=0  n 


Table  A. 7 gives  a comparison  of  the  use  of  equation  (A. 18)  and  the 
simulation  method  discussed  earlier  in  the  subsection.  The  parameter 
NK  represents  the  number  of  Monte  Carlo  runs  required  to  achieve  the 
level  of  accuracy  shown  in  the  table.  As  one  can  see,  equation  (A. 18) 
gives  the  required  solution  and  does  not  necessitate  the  substantial  run 
times  of  the  simulation. 

A. 3 LOAD  REDUCTION  PROCESS* 

A call  that  does  not  find  a free  path  in  a network  may  seize  several 
links  In  tandem  before  finally  being  blocked.  After  blockage,  the  call 

♦This  process  was  developed  by  Mr.  Dave  Garbin,  when  he  worked  at  DCEC. 
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TABLE  A. 7 LOSS  PROBABILITY  COMPARISONS  FOR  UNRELIABLE  CHANNELS  (s=20,p=10) 


.03113 

.00969 

.00511 

.00222 

.03176 

. 00946 

.00494 

.00227 

.03150 

.00963 

.00483 

.00230 

.03133 

.00964 

.00479 

.00229 

.03100 

.00963 

.00472 

.00228 

. 03098 

.00951 

.00470 

.00227 

.03082 

.00952 

.00466 

.00228 

.03067 

. 00964 

.00463 

.00228 

.03047 

.00966 

.00465 

.00229 

.03049 

.00965 

.00466 

.00228 

.03048 

.00965 

.00466 

.00228 

.03052 

.00964 

.00463 

.00228 

. 03037 

.00961 

.00460 

.00227 

. 03037 

.00959 

.00465 

.00227 

.03053 

.00958 

.00465 

.00227 

.03059 

.00960 

.00464 

.00227 

.03058 

.00961 

.00464 

.00227 

.03051 

.00962 

.00463 

.00227 

.03050 

.00958 

.00463 

.00227 

.03042 

.00957 

.00462 

.00227 

1000 

2000 

3000 

4000 

5000 

6000 

7000 

8000 

9000 

10000 

11000 

12000 

13000 

14000 

15000 

16000 

17000 

18000 

19000 

20000 


will  be  dropped  from  the  preceding  links.  Because  it  was  assumed  that 
the  set  up  time  is  zero  a blocked  call  will  not  appear  as  offered  traffic 
to  those  links  where  it  found  a free  channel,  and  the  offered  load  to 
each  link  must  be  reduced  to  reflect  this  fact. 

An  example  network  is  shown  in  Figure  A.l.  For  now,  assume  that 
there  is  one  class  of  calls  in  the  network.  The  distribution  over  all 
links  of  traffic  originating  at  the  source  node  S destined  for  node 
D.psq,  is  desired.  According  to  the  routing  table  shown  in  Figure  A.l, 
all  traffic  from  node  S passing  through  node  1 first  attempts  to  proceed 
from  node  1 directly  to  node  D.  If  any  traffic  is  blocked,  it  is  next 
routed  via  node  2 on  to  node  D.  The  traffic  initially  offered  to  each 
link  may  be  expressed  as: 


srpsD’ 

p1D=pSD^"PBS1^ 

(A. 19) 

1 2=pl DPB1 D*  and 

p2d=p12^'PB12^ 

where  PB^y  is  the  blocking  probability  on  link  XY. 
The  traffic  reaching  the  destination  node  is 


pD=p1D(1_PB1D)+p2D(1' 

pb20). 

(A. 20) 

Assuming 

the 

following  initial  values  of  link 

blocking. 

pB$i  = -05»  PB1D  = . 

1 

(A. 21 ) 

PBj2  = .15,  and  PB2D 

= .2 

and  that 

the 

offered  traffic  pcri  = 10,  it  fol 

lows  from  (A. 19)  and 

(A. 20) 

that: 

oU 

P£1=10.00,  p-j  50, 

(A. 12) 

p12=.95,  pd  =9.20  and 

p2d=.81 . 

The  traffic  1 

blocked  at  each  node  is  given  by: 

pS*pSDPBS1=’5 

(A. 23) 

°2"p2DPB2D**16 


prp12PB12sJ4  * 
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The  traffic  lost  at  node  2 (in  going  from  2 to  D);  f>2  , had  been 
carried  by  the  link  between  nodes  1 and  2.  Hence  the  initial  estimate 
of  the  offered  traffic  to  link  1,2  has  to  be  reduced  as  follows; 


p12=p12  ' • 


(A. 24) 


This  is  the  original  estimate  of  the  offered  load  on  link  1,2  minus  the 
offered  load  corresponding  to  the  traffic  that  was  blocked  on  link  2.D. 


Similarly,  for  the  link  between  nodes  S and  1,  we  have: 

P1+P2 

PS1=PS1  ■ l-PBr1  * 


(A. 25) 


Substituting  values  from  above  yields: 


p i o= • 76 


pjl  =9. 68> 


(A. 26) 


No  load  reduction  is  made  to  the  traffic  offered  the  link  between  nodes 
1 and  D as  well  as  to  the  traffic  offered  the  link  between  nodes  2 and 
0.  Neither  link  carries  any  traffic  which  eventually  is  blocked  within 
the  network;  thus, 


p 1 D~p 1 D P2D=P2D 

and 


(A. 27) 


This  procedure  is  repeated  until  the  traffic  between  all  other  source 
to  destination  nodes  has  been  distributed  (including  load  reduction), 
and  accumulated  for  each  link.  Link  blocking  probabilities  are  recom- 
puted using  the  new  estimates  of  traffic  offered  each  link  at  the  end  of 
the  load  assignment  procedure. 

For  the  case  where  one  class  of  traffic  can  preempt  another,  the 
preempted  traffic  is  load  reduced  as  if  it  were  a blocked  traffic. 
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